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ABSTRACT 

We show individual high resolution spectra of components A, B, and C of the nearby late-M type multiple system LHS 1070. 
Component A is a mid-M star, B and C are known to have masses at the threshold to brown dwarfs. From our spectra we measure 
rotation velocities and the mean magnetic field for all three components individually. We find magnetic flux on the order of several 
kilo-Gauss in all components. The rotation velocities of the two late-M objects B and C are similar (v sin / = 16km s~'), the earlier 
A component is spinning only at about half that rate. This suggests weakening of net rotational braking at late-M spectral type, and 
that the lack of slowly rotating late-M and L dwarfs is real. Furthermore, we found that magnetic flux in the B component is about 
twice as strong as in component C at similar rotation rate. This indicates that rotational braking is not proportional to magnetic field 
strength in fully convective objects, and that a different field topology is the reason for the weak braking in low mass objects. 
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1. Introduction 

o\ : 

CO . The rotation of Sun-like stars is braked fol lowing an empiri- 
,— ' ' cally determin ed braking law with v oc r 5 ( Skum anichL 1 1 972t 
<N : iBarnesLl2rjQ7l) . At several Gyr they have lost most of their an- 
\Q • gular momentum and become slow rotators like the Sun. Fully 
CD convective stars, however, are apparently not braked that much, 
C , and field dwarfs of spectral typ e late-M or L rotate more rapidly 
2 ■ than their higher mass siblings (Mohanty & BasrilPiOOll) . 
, Stars are believed to lose angular momentum due to a mag- 
■ netic stellar wind after the first phase of star formation. Angular 
r> | momentum loss critically depends on the star ' s magnetic field 
^ I and its geometry dMesj^ll984l:|j^walerlll988l:IChabover et all 
ll995HSills et al.ll2000t) . The mechanism of magnetic field gener- 
ation is of crucial importance since it determines magnetic field 
strength, hence braking, and how it reacts on angular velocity 
changes. Since the mechanism of magnetic field generation is 
believed to change at the threshold to fully convective stars, it 
can be expected that a change in the braking law appears as well. 

The lack of slowly rotating objects of spectral type later than 
mid-M could be a consequence of a change in the net braking. 
These objects, however, are so faint that it is difficult to probe 
a sample that is unbiased to luminosity effects, i.e. a luminosity 
limited sample will always contain more bright (hence young 
and less braked) low mass objects than fainter ones. 

The close hierarchical system LHS 1070 (GJ 2005) is an 
ideal probe for rotational evolution of late-M objects. Its indi- 
vidual components can be spatially resolved; the A component 
is a mid-M star (M5.5) while the two fainter components B and 
C are cooler (around M9). Measuring the astrometric orbit of 
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B and C, iLeinert et al.l (|2000) determined masses at the limit to 
brown dwarfs. 

In this paper, we measure rotation velocities and magnetic 
flux of LHS 1070 A, B, and C individually. Because all com- 
ponents are probably coeval, this gives direct insight in spectral 
type dependent braking and magnetic field generation. 

2. Data 

In 2006, the three components LHS 1070 A, B, and C were 
nearly exactly lined up on the sky, a most suitable configuration 
for simultaneous long-slit spectroscopy. We observed the triple 
in the night of October 09, 2006 during a science verification 
run (Prog. ID 60.A-9078) of the newly com missioned AO-fed 
high resolution NIR spectrograph CRIRES dKaufi et all l2006h 
at Antu (UT1) of the VLT on Cerro Paranal. 

Separations of LHS 1070AB and AC are about 1.35"and 
1.75", respectively. The spectra are spatially well resolved with 
only minimal overlap of the B and C components. Observations 
were carried out in order 57/0/i with on-target AO, total expo- 
sure time of was 480s. We chose a slitwidth of 0.4", i.e. a nomi- 
nal resolving power of 50 000. Data reduction followed the stan- 
dard procedure of dark subtraction and flatfielding. After proper 
background subtraction to remove the OH airglow emission the 
spectra where ext racted using an ad aption of the optimal extrac- 
tion algorithm of iRobertsonl d 19861) . The chosen spectral win- 
dow is clear of telluric lines on a 2% level. Thus we do not need 
to use a standard star for correction. Instead we determined the 
instrument response from the flatfield spectra and removed the 
remaining tilt by minimal rectification of the pseudo-continuum. 

Wavelength calibration was performed using a ThAr spec- 
trum in comparison to an updated line list (Kerber, priv.comm.). 
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Fig. 1. A part of our CRIRES data in the three components LHS 1070 A, B, and C (from top to bottom), smoothed by a 3 pixel 
box-car. 



Flux contamination of the B and C components was found to be 
much smaller than the intrinsic S/N ratio. Thus we do not see any 
reason to doubt on the spectral purity of any of the components 
in LHS 1070 for the subsequent analysis. A part of our data is 
shown in Fig.[T| at 9920 A the spectrum falls on a gap between 
two of the four chips. 

3. Analysis and Results 

3.1. Method 

Our observations were carried out in the Wing-Ford band 
of molecular FeH (starting at 9900 A), which shows sharp 
absorption features in mid- to late-M dwarfs. Some of the 
lines are strongly magn etically sensitive while some are not 
dReiners & Basril l2006bl) . This makes the region ideal for the 
analysis of rotational broadening and the measurement of mag- 
netic fields in ultra-cool objects. Ou r an alysis follows t he strat - 
egy used bv lReiners & Basril d2007al) and lReiners et alj d2007bl) , 

We use spectra of two template M-stars, one without any 
signs of Zeeman broadening, and the other showing strong ef- 
fects of Zeeman broadening. First, the template spectra are 
scaled according to an optical depth scaling law to fit the target's 
absorption strength. Then, the two template spectra are artifi- 
cially broadened to match the targets rotation velocity. Finally, 
we search for the linear combination of the two template spectra 
that simultaneously fits the magnetically sensitive and the mag- 
netically insensitive lines. From the interpolation parameter, we 
determine Bf, the product of the mean magnetic field and the 
filling factor, assum ing a linear relation between Z eeman broad- 
ening and Bf (see lReiners & Basril l20065l2007al) . 

As template spectra, we use data of GJ 1002 (M5.5) and 
Gl 873 (M3.5) taken with HIRES at Keck observatorjfl both 
stars have v sin i < 3 km s _1 . Because the template spectra were 
taken at lower spectral resolving power (R = 3 1 000), we applied 
a Gaussian broadening function to our data in order to match the 
resolution of the templates. This is a crucial step since it sets the 
zero-point for the calibration of rotational broadening. In order 
to check our results, we carried out the analysis with a different 
set of templates observed at higher resolution. We used a spec- 
trum of GJ 1227 (M4.5, v sin i < 3 km s _1 ) for the field-free tem- 
plate, and a sunspot spectrum for the spectrum in the presence 



of a magnetic fiekfl The spectrum of GJ 1227 was also taken 
with HIRES at Keck but at higher resolving power (R ~ 80 000). 
The template spectra were degenerated in order to match the res- 
olution of our CRIRES data. We find that both strategies provide 
highly consistent results on the values of v sin i for all three com- 
ponents of LHS 1070, i.e. differences are within lkms" 1 . For 
the following, we rely on the low resolution set of templates be- 
cause they provide a much more consistent data set, and their 
resolution is entirely sufficient for our purposes. 

3.2. Rotation 

We show our best fits over the 30 A- wide region used for fit- 
ting as red lines in all three components of LHS 1070 in Fig. [2] 
Smoothing our spectra to match the lower resolution of the 
template spectra strongly reduced the effective noise level. S/N 
ratios per resolution element are now on the order of 70 for 
LHS 1070 A, 35 and 30 for LHS 1070 B and C, respectively. 
Differences in projected rotation velocity can be distinguished 
by eye for example at 9939 A where two neighbored lines are 
just separated in the A-component but already blended in com- 
ponents B and C. The two latter objects do not show any obvious 
difference in the width of individual lines. 

We derive projected rotation velocities of v sin/ = 8, 16, 
and 16kms _1 in LHS 1070 A, B, and C, respectively. The over- 
all fit quality is rather good particularly in the A-component of 
the system. In the two late-M components B and C, the fit has 
some weaknesses which we attribute to the much stronger ab- 
sorption that may not be perfectly captured by the optical depth 
scaling. However, for our analysis of v sin i the quality is en- 
tirely sufficient. A formal x 1 analysis shows that the uncertainty 
in v sin / is on the order of 1 kms 1 for the three objects. The 
absolute uncertainty is somewhat higher since the zero-point de- 
pends on the resolution we adopt for the template spectra and 
our target spectra. We estimate that absolute uncertainties are on 
the order of 3 km s _1 , but that does not affect the comparison of 
rotation velocit i es bet ween the three components of LHS 1070. 
iBasri & Marcvl d 19951) have measured v sin / in a spatially unre- 
solved spectrum of LHS 1070. They find v sin / = 8 + 2kms _1 
in excellent agreement with our result for LHS 1070 A, which 
dominates the spatially unresolved spectrum. 



1 We thank G. Basri for providing the template spectra. 
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Fig. 2. Data (black) and fit (red) in the entire region we used for 
the fit in the three components LHS 1070, A, B, and C (from top 
to bottom). The data was degenerated in resolution to match the 
resolution of the template spectra (R — 31 000). 



3.3. Magnetic flux 

Magnetic flux was determined by searching for a linear inter- 
polation between a field free template spectrum and a template 
spectrum affected by magnetic Zeeman broadening. For the lat- 
ter, we use a spectru m of Gl 873 and we ad opt a magnetic flux 
of Bf = 3.9 kG (see lReiners & Basrill2007al) . This implies that 
our results cannot exceed values of ~ 4kG. 

At every wavelength, the interpolation lies between the two 
template spectra. In magnetically insensitive regions, the dif- 
ference is very small while at wavelengths that are affected by 
Zeeman broadening, a difference proportional to the magnetic 
flux occurs. We show a representative region of our spectra at 
the wavelength range 9946-9956 A in Fig. [3] Here, we highlight 
three absorption lines that are particularly sensitive to Zeeman 
broadening. We plot the target spectra in black and overplot the 
(rotationally broadened) template spectra of GJ 1002 (no mag- 
netic flux) in blue, and the one of Gl 873 (strong flux) in red. 
The interpolation that best fits the data is overplotted in green. 

A very accurate fit is achieved in the spectrum of 
LHS 1070 A, from which we derive a magnetic flux of 2.0kG. 
The fit quality in the spectrum of LHS 1070B is not as ac- 
curate. However, it appears that the two sensitive absorption 
lines at 9948 A and 9954 A are strongly smeared out, possi- 
bly even more than in the template spectrum of Gl 873. This 
indicates very strong magnetic flux, and we adopt a value of 
Bf ~ 4 kG for component B (note that the red line is covered 
by the green fit in this panel). The presence of strong magnetic 
flux can also be seen by comparing the spectrum of LHS 1070 B 
to the one of LHS 1070C in the third panel. Here, the fit qual- 
ity is much higher again, we deduce a magnetic flux of 2kG 
as in the A-component. Comparing this spectrum to the one of 
LHS 1070B (particularly the lines at 9948 Aand 9954 A), it ap- 
pears that the magnetically sensitive absorption lines are much 
weaker in LHS 1070 B than they are in LHS 1070 C 

As discussed in iReiners & Basril d2006bl l2007al) . our lack 
of knowledge about the Zeeman splitting patterns of the FeH 
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Fig. 3. Portion of the spectral range that shows the effect of mag- 
netic flux. Green bars mark magnetically most sensitive lines. 
Top to bottom: LHS 1070 A, B, and C. Black lines are the target 
spectra, red lines mark spectra affected by strong magnetic flux 
(G1873), blue lines spectra without magnetic flux (GJ 1002). 
Green lines are our best fits. Note that in the B component the 
red line is covered by the green one. 

Table 1. Parameter of the t hree components LH S 1070 A, B, and 
C. Spectral types are from lLeinert et alj d2000t) . 



Component 


SpType 


V 


J 


v sin z 


Bf 










[km/s- 1 ] 


[kG] 


LHS 1070 A 


M5.5 


15.42 


9.25 


8 


2.0 


LHS 1070 B 


M8.5 






16 


4.0 


LHS 1070 C 


M9.0 






16 


2.0 



molecule causes the main uncertainty in the determination of 
magnetic flux using this molecule as a tracer. On an absolute 
scale, we estimate our uncertainties in magnetic flux on the order 
of a kilo-Gauss. However, the main result of our magnetic flux 
analysis is already seen in a rough comparison of the three spec- 
tra: All three components show a substantial effect of magnetic 
flux on the spectral lines, and the magnetic flux in LHS 1070B 
is stronger than the flux in LHS 1070 A and LHS 1070 C. 

4. Summary and Conclusions 

We have isolated high-resolution spectra of the three compo- 
nents of the multiple M-object system LHS 1070. From each 
spectrum, we have determined the projected rotation velocity 
v sin and the mean magnetic field over the entire star. The results 
of our analysis are summarized in TableQ] We find a projected 
rotational velocity of v sin i = 8 km s _1 in the hottest component, 
LHS 1070 A, while rotation in LHS 1070 B and C is about twice 
as rapid. Under the reasonable assumption that the entire system 
LHS 1070 was formed at the same time, we are facing three dif- 
ferent objects that have evolved from the sa me initial conditions . 
Disc orientations in pre-main sequence stars (Monin et al., 2006) 
and orbit orientations in multiple systems (ISterzik & TokovirrrnL 
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Fig. 4. Projected rotation velocities of LHS 1070 A, B, and 
C. We overplot rotation velocities according to a modified 
Skumanich braking law at 1, 2, 5, and and lOGyr (top to bot- 
tom). 



120021) are partially correlated, and iBate et al.l (l2000h find that 
strong misalignments are unlikely in binaries with separations 
< 100 AU. Thus, we may assume that the inclin ation angles i 
of LH S 1070 A, B and C are comparable (see also lLeinert et all 
1200 11) . As long as they suffer the same rotational braking, we 
would expect comparable rotation velocities in all components. 
We observe comparable rotation rates in the components B and 
C, which are of very similar spectral type. However, the rotation 
velocity in LHS 1070 A is about a factor of two smaller. This 
can probably be attributed to a difference in the rotational brak- 
ing du ring their evolu tion. 

iMohantv & B~asrH d2003l) showed that in the field the frac- 
tion of slow rotators (v sin; < 3kms"') later than spectral type 
M6 is much smaller t han among early-M d warfs, which gener- 
ally rotatej[ery_slowly (Delfo sse et all [T998). However, the sam- 
ple of IMohantv & Basril d2003l) could be biased to bright (hence 
young) objects that are still more rapidly rota t ing tha n their 
older (hence fainter) siblings. iReiners & B asri (2006a) found 
slow rotation in a mid/late-M subdwarf, and very rapid rotation 
(v sin; = 65kms~ I ) in a late-L subdwarf. Since late-type sub- 
dwarfs are probably the oldest relics of the early galaxy having 
spun down for several Gyr, this is another argument for less effi- 
cient rotational braking at later spectral type. 

In LHS 1070, we see three objects of different spectral type 
that are probably of same age. We plot the projected rotation 
velocities of LHS 1070 A, B, and C in Fig.|4] As stated above, 
the two later objects show higher rotation rates than the earlier 
A-component. One explanation may be a difference in th e ini- 
tial conditions, but together with Mo hantv & Basril (120031) this 
would suggest that later objects in general start off rotating more 
rapidly, which is unlikely. On the other hand, the rotation activ- 
ity connection seems to vanish around late-M spectral type. This 
is likely connected to a modified wind-law and consequently to 
different rotational braking. 

Detailed simul a tions o f rotational brakin g were performed 
by IChaboveretaTJ (119951) : ISills et al l d2000t) . They give a pre- 
scription for a wind-braking law in stars at given angular veloc- 
ity, and they show that the braking law in fully convective stars 
requires different treatment (to fit ob servations in the Hyades 
they had to tune their value of w cr i t ; ISills et all 120001) . In this 
work, we do not intend to give a full description of the braking 
law in fully convective objects. Instead, we assume that the net 



braking law in LH S 1070 A is still Skumanich-type (v oc f -0 - 5 , 
Skumanichl |1972|) . and that the exponent of the braking law is 
gradually decreasing towards lower temperatures. 

For this purpose, we calculate the rotational evolution start- 
ing at an age of lOMyr and an initial rotation velocity of 
70kms~'. For our spectral type dependent braking law, we chose 
v oc r a with gradually decreasing a = 0.5-0.3 in the spectral 
range M5-M9. Our results are overplotted as dashed lines in 
Fig-El Four lines for ages of 1, 2, 5, and lOGyr are shown (top 
to bottom). The assumed spectral type dependent scaling law 
can reproduce the difference in rotation velocities as measured 
in LHS 1070, and at the same time generate the rise in mini- 
mum rotational velocities con sistent to measurements in late-M 
stars (Mohantv & Basri, 2003). We thus conclude that braking 
in late-M stars is spectral-type dependent, and that the lack of 
slowly rotating late-M dwarfs is not a selection effect due to a 
bias towards young rapid rotators. 

The age estimated from our approach is on the order of 1 Gyr. 
iBasri&Marcvld 1995b found a high velocity component perpen- 
dicular to the galactic plane that suggests higher age. We think 
that this point deserves further attention. Although age estimates 
from space velocity are rather uncertain, this may hint to an ab- 
solute scale of the braking law different to the one used here. 

All three components of LHS 1070 show strong magnetic 
flux suggesting that in fully convecti ve objects, the field s trength 
does not strongly depend on rotation. fLeinert et al.l (12000) found 
signs of magnetic activity in components A and B, but not in 
C. This is consistent with our results in the sense that the mag- 
netic flux in LHS 1070B is much stronger than the one in 
LHS 1070 C. Nevertheless, one expects that the C component 
also exhibit s Ha emission, whic h may still be below the detec- 
tion limit of iLeinert etail (120001) . 

Interestingly, the difference in magnetic flux between 
LHS 1070 B and C does apparently not influence rotational brak- 
ing. It is generally believed that angular momentum loss de- 
pends on magnetic field s t rength with an exponent given by 
the field geometry (Mestel, 1984). This c oupling, howeve r, de- 
pends on the field generation mechanism (iKawaleij 1988), and 
angular momen t um lo ss becomes s aturated at a certain leve l 
(IChaboveretallll995h . According to lChabrier & Kiikerl (l2006). 
large scale magnetic fields of equipartition strength (a few kG) 
can be generated by an a 2 dynamo. These authors predict the 
field topology to differ from an organized dipole field, which 
could explain our result that rotational braking does not grow 
with magnetic field strength. It would imply that field geometry 
is the reason for lower rotational braking in low mass objects. 
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